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Abstract 

Hydrogen reduction of the lunar regolith simulants JSC-1A and LHT-2M is investigated in this paper. 
Experiments conducted at NASA Johnson Space Center are described and are analyzed utilizing a 
previously validated model developed by the authors at NASA Glenn Research Center. The effects of 
regolith sintering and clumping, likely in actual production operations, on the oxygen production rate are 
studied. Interpretations of the obtained results on the basis of the validated model are provided and linked 
to increase in the effective particle size and reduction in the intra-particle species diffusion rates. Initial 
results on the pressure dependence of the oxygen production rate are also presented and discussed. 

Nomenclature 

c a molar concentration of hydrogen in regolith bed (moles/m 3 ) 

D effective gas diffusion coefficient (m 2 /s) 

F factor defined by Equation (1) (s -1 ) 
k equilibrium constant 
L length of regolith bed ( m) 
p gas pressure (Pa) 
r p regolith particle radius (m) 
t time (s) 

t c time for complete conversion (s) 

T temperature (K) 
u velocity through reactor (m/s) 
a particle conversion fraction 

p a molar concentration of iron oxide in regolith (mol/m 3 ) 
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I. Introduction 


A major objective of lunar in-situ resource utilization (1SRU) (Refs. 1 to 3) is production of oxygen 
from lunar regolith. Mare basalts containing oxygen include olivine, ilmenite, and pyroxene (Ref. 4). 
Energy has to be supplied to remove the oxygen, i.e., the associated chemical reactions are endothermic. 
Free energy considerations show that the iron containing minerals require the least energy input for this 
purpose (Ref. 5). Various technologies have been proposed for the chemical conversion. These include 
ilmenite reduction by hydrogen, molten electrolysis, and carbothermal reduction using methane (Refs. 6 
and 7). The hydrogen reduction method targets the iron oxide in the regolith. Iron oxide occurs in lunar 
regolith as ilmenite, particularly in titanium-rich basalts, or may be found dispersed in the glassy phase of 
lunar regolith particles. The overall gas-solid reaction between FeO and hydrogen is given by: 

FeO + H 2 <=> Fe + H 2 0 (A) 

The formed water may be electrolyzed in an electrolyzer subsystem to produce oxygen and recycle the 
hydrogen. 

Hydrogen reduction systems are well developed for Earth applications but mainly for the purpose of 
metal extraction, as opposed to water vapor and oxygen extraction for lunar applications. In addition, in 
the lunar context, the reactor must interface with the other sub-system processes such as upstream regolith 
extraction and beneficiation and downstream electrolysis and phase separation in a way that establishes 
the most favorable balance between efficiency, robustness, and equivalent system mass. A need exists, 
therefore, to establish an analytical framework to couple these sub-systems together. The building blocks 
of this framework are the individual sub-system models. The development of a generalized analytical 
model of the chemical reactor sub-system and its validation against experimental data from reduction of 
JSC-1A, a lunar regolith simulant, has been discussed in an earlier paper (Ref. 8). 

The primary objective of this paper is to investigate the hydrogen reduction of a new simulant, 
LHT-2M, which has been manufactured to resemble the highland type of lunar regolith, and to compare it 
with the hydrogen reduction of JSC-1A. To this end, experiments have been conducted at NASA Johnson 
Space Center for analysis with the validated model developed at NASA Glenn Research Center. This 
model (Ref. 8) is formulated in terms of the following overall steps: (i) inflow of pure hydrogen into the 
regolith containing reactor, (ii) reaction of hydrogen with the iron oxide portion of the regolith, and (iii) 
outflow of left-over hydrogen and produced water vapor. Step (ii) is the gas-solid reaction step which is 
the key to the conversion. Results of regolith particulate analysis conducted on Apollo-returned samples 
(Ref. 6) suggest that the particles are not dense solids but have some porosity that would enable gases to 
diffuse into the particles, i.e., the iron oxide portion can be essentially completely converted with the 
reduced iron remaining behind in the particulates. It may also be noted that both JSC-1A and LHT-2M 
have significant glass content with iron oxide dispersed in the glass phase. The diffusion of gases in the 
glass phase is expected to be the controlling feature for the reduction process. These observations support 
utilization of a “shrinking core” diffusion-reaction model (Ref. 9) for the gas-solid reduction process for 
the simulants. This model should be valid for lunar regolith as well if the gas diffusion in the particles 
remains the controlling step. The reaction model is coupled with the inflow/outflow of gases to complete 
the reactor sub-system model. 

Two key parameters appear in the model: (i) a diffusion time related to the diffusion of gases in the 
interior of the regolith particles, and (ii) a residence time of the gases in the reactor. It is found that for 
high levels of conversion or for a relatively small ratio of residence time to diffusion time, the governing 
parameter is the diffusion time. Validation of the model was conducted by comparison to experiments 
carried out at Pacific Northwest National Laboratory (PNNL) (Ref. 10) and at Lockheed-Martin (Ref. 11) 
where the quantity of processed regolith simulant in the experiments was about three orders of magnitude 
larger than that used at PNNL. 
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A factor F, which is inversely related to the diffusion times for different particle sizes, is derived from 
the model formulation (Ref. 8) and is defined as follows: 


F = 6— 
Pa 


k D 


( 1 ) 


It was found, using the PNNL data for JSC-1A in conjunction with the model, that the factor F may be 
correlated to the temperature by an Arrhenius expression: 

log(F) = - 10 ^ 61 - 2.8808 (2) 

The above expression appears to work well for JSC-1A in the temperature range of 873 to 1273 K, i.e., in 
between 600 and 1000 °C, under atmospheric pressure conditions. 

The factor A depends upon the concentration of hydrogen, c„, in the reactor, the intra-particle species 
diffusivity, D, the equilibrium constant, k, for the iron oxide reduction reaction, the concentration of iron oxide 
in the regolith/simulant, p„, and the effective particle radius, r p . The concentration, c 0 , of hydrogen, provided 
that it is the sole inflow species, is proportional to the reactor pressure,/!, assuming perfect gas behavior. The 
diffusion coefficient, D, will depend upon the regolith/simulant of interest and is, in general, a function of both 
reactor temperature, T, and pressure,/!. The equilibrium constant, k, has generally been presumed to be a 
function of temperature alone and to be independent of pressure for the reduction of iron oxide with hydrogen 
(Ref. 8). Both p a and r p will depend upon the characteristics of the utilized regolith/simulant. It is clear 
therefore, that much work remains to be done to further characterize the variations in the factor F for different 
regoliths and simulants and to understand the impact of reactor pressure. 

It is also found that, apart from the dependencies noted above, regolith/simulant sintering and 
clumping during processing can impact the water production rate. A quantification of this effect is 
provided. In addition, some results on the pressure dependence of F, and hence of production rate, are 
also presented and discussed. In the following, a brief discussion of the developed model is provided first, 
and a description of the experimental apparatus and procedure is given next. This is followed by a 
presentation and discussion of the results. 


II. Model Description 


The shrinking core formalism (Ref. 9) is utilized in the model and is schematically depicted in 
Figure 1. Conversion proceeds from the exterior particle surface inward. As the reaction progresses into 
the interior of the particle, the reactant and product gases must diffuse through the particle matrix 
structure formed after the conversion. 



Unreacted core 


Particle radius, !/ 

Figure 1. — Shrinking core model configuration. 


NASA/TM- 


- 2011-216993 


3 



The equation for the conversion fraction a (related to the fraction of iron oxide converted) at a given 
time, t, is given by (Ref. 8) 


3[l-a] 2/3 -2[ 1 -a]- 1 


Ft 


1 + 3(l - s) 

( \ 

D 

2 


VP J 


[I -a r 

1 - (l - a) 1 


( 3 ) 


where the factor F is given by Equation (1), L is the length of the reactor bed, and u is the gas velocity 
through the reactor. For the case when the residence time, Liu, of the gases is small compared to the 
diffusion time, r p !D, of gases in the particle matrix the equation for a becomes (Ref. 8) 

3[l-a] 2/3 -2[l-a]-l = -Ft (4) 

It may be noted that the time required for high conversion (i.e., as a approaches unity) is also given 
by Equation (4). In practice, this means that unless the residence time is large compared to the diffusion 
time, the time for complete conversion is governed only by the diffusion time of the gases inside the 
particles. This implies that the factor F is a characteristic descriptor of the reduction process. Examination 
of the expression for F (see Eq. (1)) allows for the following observations: 

(a) F has units of inverse time (i.e., s _1 ). 

(b) the equilibrium constant, k, is for Reaction (A) which is, presumably, not simulant dependent. 

(c) the FeO concentration, p„, is simulant dependent but should be relatively invariant from sample to 
sample or relative to the particle size of the same simulant. 

(d) the hydrogen concentration, c 0 , is relatively uniform in the reactor when the residence time is small 
compared to the diffusion time and is an operator controlled parameter of the experiment; it is a function 
of the reactor pressure, p, and temperature, T. 

(e) the particle size, r p , is an effective particle radius. There is generally no ambiguity about particle size 
for a relatively narrow size range; however, when the particle size range is broad, experiments need to be 
done to characterize a physically meaningful “effective” particle size. 

(f) if the factors in (b) to (e) above are known, then F may be directly related to D, the effective diffusion 
coefficient of the gases in the particle by Equation (1). 

III. Experiment Description 

An image of the reactor schematic and experimental set-up are shown in Figure 2. The reactor tube 
has an internal diameter of 2.12 cm, an overall length of 76.2 cm, and a test section (filter-to-filter) length 
of 22.9 cm. It is constructed of Inconel 625 and is enclosed using flared T-fittings. A 50 pm ceramically 
bonded alumina disk (2.08 cm outer diameter by 0.32 cm inner diameter by 0.635 cm thick) is attached to 
a 0.32 cm thermocouple using compression ferrules to keep it from moving. The thermocouple and 
attached filter are inserted through the bottom T-fitting. A similar sintered metal filter is attached to a 
0.32 cm rod through the top flared T-fitting. Together they serve to keep particulates from flowing out of 
the reactor and hold the simulant bed in place while allowing gas to pass through. Gas flows into and out 
of the reactor using the 90° legs of the top and bottom T-fitting, respectively. The simulant is heated 
through the reactor walls using an external tube furnace. Temperature control is achieved using the 
internal thermocouple with the attached filter. The reactor is designed to withstand a temperature of 
1323 K (1050 °C) at an absolute pressure of 4.4 atm. Two mass flow controllers are used to set the 
desired gas flow rates of inert gas and hydrogen. Pressure transducers are installed on the inlet and outlet 
sides of the reactor to measure the pressure and pressure drop in the reactor. A flow meter and a back 
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Figure 2. — Reactor schematic and experimental set-up. 


pressure regulator placed downstream of the furnace complete the flow system. A relative humidity 
sensor, with an accuracy of 1 percent RH, is used to monitor the water concentration in the outflow gas. 
Water is removed from the outflow gas in a dryer. All data are recorded digitally. 

A batch size of approximately 10 g of simulant is used for each test. A gas flow velocity of 1 cm/s (as 
measured under atmospheric conditions) is utilized. Previous experiments suggested that this flow rate is 
generally sufficient to keep the simulant agitated and fluidized which is useful in preventing 
agglomeration and sintering of the particles during the test. However, sintering and clumping of the 
particles may occur if the fluidization is not isotropic, e.g., due to flow non-uniformities. The 
temperatures at which the reduction is carried out ranged between approximately 1023 to 1323 K (i.e., 

750 to 1050 °C). There is usually some bound water in the simulant sample which must be baked off prior 
to reduction in order to avoid ambiguity between reaction produced water and bound water. Thus, the test 
procedure is as follows: 


1 . Load reactor with simulant. 

2. Flow inert gas, typically helium, while raising simulant to the test temperature. 

3. Switch from inert gas to hydrogen flow after simulant has reached test temperature and evaporation of 
bound water is complete as evidenced by the humidity sensor data. 

4. Continue hydrogen flow until weight loss from water production during reduction is below the 
humidity sensor accuracy threshold or for a maximum of 3 hr whichever comes first. 

5. Assess sintering and/or clumping of particles after test is complete. 

The water output as a function of time is derived from the output of the relative humidity sensor. The 
dryer is weighed before and after heating the simulant to obtain the bakeout water amount. The dryer is 
also weighed after the completion of each test to obtain the overall amount of water produced by 
hydrogen reduction during the test. 
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IV. Results and Discussion 


A. Chemical Composition and Size Distribution Comparison of JSC-1A and LHT-2M 

The major constituents of JSC-1A and LHT-2M by weight percent, as obtained from their material 
safety data sheets, are shown in Table 1 and can be seen to be similar. A notable difference is the amount of 
iron oxide in the two simulants which is of interest in this study since hydrogen reduces mainly iron oxide in 
the temperature range of interest. The iron oxide generally consists of both FeO and Fe 2 C> 3 . The maximum 
amount of water that can be obtained by hydrogen reduction of the oxides may be estimated based upon this 
composition. This is found to be approximately 2.8 and 1.2 wt% for JSC-1A and LHT-2M, respectively. It 
may be noted that this maximum amount of water that can be obtained from the hydrogen reduction of iron 
oxide corresponds to a conversion fraction, a, of unity in Equations (3) and (4). 

The particle size distributions of JSC-1A and LHT-2M, obtained from previous size sorting studies at 
NASA, are shown in Figure 3 which plots the cumulative mass percentage as a function of particle size. It 
is seen that the size distribution of the two simulants are very similar. 

The similarities in the major chemical constituents of JSC-1A and LHT-2M suggests that the intra- 
particle species diffusivity may also be similar since hydrogen would need to diffuse through similar 
compositions (and, presumably, internal structures) to reach the dispersed iron oxide. In addition, the 
particle size distribution similarity suggests that the effective particle size, r p , of both simulants should 
also be similar. The quantity, p a , in Equation (1) is approximately a factor of three smaller for LHT-2M 
according to the reported concentration of iron oxide (Table 1). Referring to Equation (1), it is then 
reasonable to expect that under the same conditions of reactor temperature and pressure, the quantity F 
will be of the same order for both JSC-1 A and LHT-2M. 


TABLE 1.— CHEMICAL COMPOSITION OF JSC-1 A AND LHT-2M 


Composition 

Weight percent JSC-1A 

Weight percent LHT-2M 

Silicon dioxide 

46 to 49 

46.7 

Aluminum oxide 

14.5 to 15.5 

24.4 

Calcium oxide 

10 to 11 

13.6 

Magnesium oxide 

8.5 to 9.5 

7.9 

Iron oxide 

10 to 11.5 

4.16 

Sodium oxide 

2.5 to 3 

1.26 

Potassium oxide 

0.75 to 0.85 

0.08 

Manganese oxide 

0.15 to 0.2 

0.07 

Titanium dioxide 

1 to 2 

0.41 

Phosphorous pentoxide 

0.6 to 0.7 

0.15 



Figure 3. — Cumulative mass percentage of JSC-1A and LHT-2M as a function of particle size. 
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B. Effects of Sintering/Clumping on Reaction Rate 

In assessing the quantities that affect the factor F and the water output rate, the sintering and/or 
clumping of the simulant particles has not been quantitatively addressed to date. Studies of sintering of 
lunar regolith and simulants have shown that softening of glass particles is the primary cause of sintering 
(Ref. 12). Sintering welds particles together due to viscous relaxation of the glass in the sample. As 
temperature increases, glass viscosity and time for sintering decrease. Experiments with Apollo 14 
composition glass demonstrated (Ref. 13) that sintering can occur within minutes at temperatures of 
1023 K (800 °C) if the glass particles are otherwise undisturbed (e.g., by flow). Fluidization of the 
regolith bed decreases the amount of sintering as particles do not remain in contact with each other long 
enough for sintering to occur. It has also been found that smaller particles tend to sinter at earlier times 
compared to larger particles. This is summarized in Figure 4 which is taken from Cooper (Ref. 12) 
following Wosinski et al. (Ref. 13) 

Referring to Equation (1), sintering maybe expected to modify both the effective particle size, r p , and 
the effective species diffusivity D. It is clear that the effective size will increase as particles stick together 
and clumps are formed. As r p increases the value of F will decrease. The impact on D is more involved. The 
sticking of particles likely alters the path through which diffusion occurs tending to reduce the “porosity” 
and making the particles more impervious to mass transfer. This would indicate a reduction in D and, 
consequently, a reduction in the value of F. Thus, the overall effect of sintering is expected to result in a 
decrease in the observed value of F, which, in turn, implies a reduction in the water production rate. 

This assessment has been verified during the current experiment series. Though not expected based 
upon the discussions presented in Section IV. A., tests with LHT-2M at temperatures above 1 173 K (i.e., 
900 °C) with batches that included particle sizes smaller than 125 pm yielded values of A that were 
substantially lower than those for correspondingly sized JSC-1 A particles. Examination of these LHT-2M 
samples after reduction showed several sintered clumps, some of which were bigger than 1 cm in size. On 
the other hand, for unsieved LHT-2M which consists of a significant fraction of larger size particles (see 
Fig. 3) and for lower temperatures, the values of F for JSC-1A and LHT-2M were comparable. Table 2 
compares the values of F for LHT-2M and JSC-1 A for several cases showing this trend. Images of three 
of the post-reduction LHT-2M samples are shown in Figure 5 for both sintered and unsintered cases. 

Temperature, °C 



5 7 9 11 13 15 17 19 

Log viscosity (poise) 


Figure 4. — Time for incipient sintering between two glass spheres deforming 
by viscous flow. (From Cooper (Ref. 12)). 
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TABLE 2.— COMPARISON OF THE FACTOR F FOR JSC-1 A AND LHT-2M FOR 


DIFFERENT PARTICLE SIZE RANGES AND TEMPERATURES 


Particle size range 

Reduction temperature, 
K 

F for JSC-1 A, 
sec -1 

F for LHT-2M, 
sec" 1 

<37 gm 

1273 to 1323 

6.7e" 5 

5.0s~' (sintered) 

40 to 75 grn 

1273 to 1323 

1.7e" 5 

9.0e"' (sintered) 

Unsieved 

1273 to 1323 

l.le" 5 

l.le" 6 (sintered) 

Unsieved 

1173 to 1200 

6.7e" 6 

2.8e" 6 (partially sintered) 

Unsieved 

1048 to 1073 

2.3e" 6 

3.1e 6 (unsintered) 



-t-es-sr sr/ 2 - 7 /'” 

, C2* — 

X 1. t'-U 5r I 



(a) (b) (c) 

Figure 5. — Post reduction simulant samples (a) <37 pm (b) 40 to 75 pm (c) 125 to 250 pm. Sintering 
and clumping is observed for (a) and (b) but not for (c). Reduction temperatures were between 1273 
and 1323 K at 1 atm reactor pressure. 


The sintering of LHT-2M in the current tests is believed to be due to non-isotropic fluidization of the 
simulant bed as was observed in other studies (Ref. 12). In order to verify this hypothesis, new 
experiments were conducted with JSC-1A which, in earlier tests using experimental hardware with 
different types of fluidization (Ref. 8), did not exhibit sintering up to temperatures of 1273 K (1000 °C). 
Two tests were conducted with JSC-1A in the present set-up at reaction temperatures of 1273 K 
(1000 °C) and 1073 K (800 °C), respectively. The post-reduced JSC-1 A sample from the 1273 Ktest 
exhibited the clumping and sintering that was observed with LHT-2M. Figure 6 shows the experimental 
conversion with time for this test with the sintering effect present. Also shown in this figure are the 
expected range of conversion from our validated model predictions and data from a previous experiment 
conducted at Lockheed-Martin. The Lockheed-Martin experiment was a rotating reactor configuration in 
which the simulant did not exhibit sintering. It is seen that for the experiment conducted during the 
current study, the conversion rate (i.e., fraction converted in a given amount of time) is about 50 percent 
smaller than expected without sintering. 

In contrast to the test case conducted at 1273 K, the post-reduced JSC-1 A sample from the 1073 K 
test did not have any observable sintering. The results from this test are shown in Figure 7. Also shown, 
as before, are the expected range of conversion from our validated model predictions as well as 
corresponding data from the Lockheed-Martin experiment. It is seen that the conversion rate from the 
current experiment agrees satisfactorily with the model estimate and prior experimental result. 

While in the past, it has been suspected that sintering would result in lower conversion, this is the first 
time that this effect has been shown quantitatively and points to the importance of maintaining good 
fluidization and gas-solid particulates mixing in the reactor. It may be noted that non-uniformities in the 
inlet flow in the current reactor are believed to cause the non-isotropic fluidization and may be caused by 
blockages in the sintered plate of the inlet section. In addition, the observations here are consistent with 
previous findings (e.g., see Fig. 4) that sintering is promoted at the higher temperatures. 
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Figure 6. — The effect of sintering on the conversion of JSC-1 A at 1273 K. 



Figure 7. — Conversion of JSC-1 A at 1073 K. Effects of sintering are minimal. 
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C. Effect of Reactor Pressure 


Reactor pressure is expected to impact the factor F through the concentration of hydrogen in the 
reactor, c a , and the intra-particle species diffusivity D as may be seen from Equation (1). In addition, it 
may be noted from Equation (4) that the time for complete conversion (i.e., a = 1) is given by 


For given simulant characteristics and reaction temperature, and noting that the equilibrium constant may 
be assumed to be independent of pressure, as noted in the Introduction section, it is clear that F is 
proportional to the product of the hydrogen concentration, c a , and the effective species diffusivity, D 

F oc c 0 D (6) 

so that the time for complete conversion becomes inversely proportional to this product 

*c K ( c o D )~ 1 (7) 

Since the inlet hydrogen concentration, c a , is proportional to the reactor pressure, the conversion time 
becomes inversely proportional to the product of the pressure, p, and the diffusivity, D 

t c cc(pD)- 1 ( 8 ) 

The effect of reactor pressure for levels greater than atmospheric has been investigated in the current 
set of experiments using LHT-2M simulant particles in the 125 to 250 pm size range. As noted earlier 
(e.g., see Fig. 5), particles in this size range did not exhibit sintering even at temperatures as high as 1273 
to 1 323 K. This was found to be true even at pressures up to 3 atm which was the maximum reactor 
pressure tested. The experimentally observed conversion rates at 1 and 3 atm pressure and 1300 K are 
plotted in Figure 8 during the initial period of conversion for the larger, unsintered particles in the size 
range of 125 to 250 pm. Until approximately 0.25 hr, the conversion from the lower pressure case is 



Time, hrs 

Figure 8. — Measured conversion as a function for 125 to 250 gm size LHT-2M at 1 and 3 atm 
pressure at a temperature of 1300 K. 
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greater than that from the higher pressure case. Thereafter, the conversion from the higher pressure case is 
greater. In fact, for the lower 1 atm pressure case, the readings from the humidity sensor go below its 
sensitivity threshold at a conversion time of approximately 0.3 hr. This behavior may be explained by 
noting that during the initial period of conversion, the iron oxide in the surface layers of the particles are 
reduced. Diffusion rates in the bulk gas phase and inside the surface pores of the particles are likely 
important during this time and these are known to be inversely dependent upon the pressure. 

The factor F is estimated from the slope of the conversion versus time curve toward the end of the 
measurement period as has been discussed previously (Ref. 8). It is found that F is about 50 percent larger 
for the 3 atm case compared to the 1 atm case. From the equations above, the pressure dependence of D 
for this case may be shown to be 


r-v —0.63 sr\\ 

D oc p (9) 

A possible explanation for this finding is obtained by considering the mechanism of species diffusion 
in the interior of the simulant particles. If the diffusion process inside the particle is governed by gas 
phase processes alone D would be inversely proportional to p. On the other hand, if diffusion of the 
species is through migration in the solid phase, the reactor pressure would not be expected to play a role 
implying that D would be pressure independent. The p 0 63 dependence of D, therefore, suggests that the 
species diffusion inside the simulant particles shares features from both gas phase and solid phase. This 
result is preliminary and more data are required at different particle sizes, reactor temperatures, and 
pressures to better establish the pressure dependence of the hydrogen reduction process. In this regard, it 
may be noted that in previous work at Lockheed- Martin, experiments with JSC-1A showed that the 
conversion rate was effectively doubled as reactor pressure increased from 0.5 to 1 atm. Thus, at sub- 
atmospheric reactor pressures, the conversion time is proportional to p which from the equation above 
suggests that the diffusivity D is largely pressure independent, presumably due to being dominated by 
solid-phase diffusion. 


V. Conclusions 

The hydrogen reduction characteristics of LHT-2M, a lunar highlands regolith type simulant has been 
studied and compared with those of the simulant JSC-1A. Experiments were conducted in a tubular 
furnace reactor. Parameters varied were particle size, reaction temperature, and reactor pressure. 
Consistent with earlier studies on lunar regolith sintering, it was found that smaller particles and higher 
temperatures promoted clumping and sintering. The sintering/clumping resulted in significant (as much as 
50 percent) reduction in the reaction rate. The deleterious effects of sintering were explained in terms of a 
model previously developed by the authors at NASA Glenn Research Center and were related to increase 
in the effective particle size and decrease in the effective intra-particle species diffusivity. The effects of 
reactor pressure on the reaction rate were also investigated. It was found that for pressures in the range of 
1 to 3 atm, the effective intra-particle species diffusivity varied as p °' 63 . The results of the present 
investigation will be useful for the design of proposed hydrogen reduction reactors. 
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